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NOZZLE FLOW WITH VIBRATIONAL NONEQUILIBRIUM 


Introduction 

We consider the problem of modeling a high temperature diatomic gas iV 2 flowing 
through a converging-diverging high expansion nozzle. The problem is modeled in two 
ways. The first model uses a single temperature with variable specific heats els functions 
of this temperature. For the second model we assume that the various degrees of freedom 
all have a Boltzmann population distribution which means that each degree of freedom 
has its own temperature and consequently each system state can be characterized by these 
temperatures. This suggests the formulation of a second model with a vibrational degree 
of freedom as having its own temperature along with a rotational-translational degree of 
freedom with its own temperature. Initially the vibrational degree of freedom is excited 
by heating the gas to a high temperature. As the high temperature gas expands through 
the nozzle throat there is a sudden drop in the rotational-translational temperature along 
with a finite relaxation time for the vibrational degree of freedom to achieve equilibrium 
with the rotational- translational degree of freedom. That is, the temperature change that 
occurs when the iV 2 gas passes through the nozzle throat is so great that the changes in 
the vibrational degree of freedom lags behind the rotational-translational energy changes. 
The resulting relaxation time is finite. It is in this context that the term nonequilibrium 
is used. That is, the term nonequilibrium denotes the fact that the energy content of the 
various degrees of freedom are characterized by two temperatures. We neglect any chemical 
reactions resulting from the high temperatures which could also add nonequilibrium effects. 

We develop the basic equations for the two models in various forms in order to check 
the derivations with other sources, references [1],[2]. The final form which is solved nu- 
merically are the scaled equations in a conservative dimensionless form. 


Single Temperature Model 


For our first model we assume that there exits a single temperature T which charac- 
terizes the energy state of the system. Using the list of symbols given in the Appendix 
A, the basic equations describing the flow through a nozzle with cylindrical symmetry are 
given by 


1 / P 

Continuity + gV • V = 0 

DV 

Momentum f>-^~ = +V • (P) 

Pij = -PSij + v(Vi,j + V jti ) + a Sij(y • V) 

de - -# dQ 

Energy g— + gV • Ve + PV • V + V • q = — + $ 
at at 

Equation of State P = gRT 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 
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where ^ -j- V ■ V is the material derivative and 


<1 — Qr&r "f~ Qz^zy 

.dT 


V = V r e r + V z e 2 
dT 


q ' = -%■ «* = az 


C„ = C^rt + C^u, Curt — 5R/2, C, 

rT c R<f) 

u v ai = - i o)+ — 

'T 0 


-tf) 




/• T 5 

;(T) = / C„dT= -R(r-T 0 ) + 
./r 0 2 


( e ^/T_i)2 

R<f> 


e^/T — i e^/^o — 1 


(6) 

(7) 

<f> — hu/k ( 8 ) 

(9) 


We assume that the external- heat sources Q are zero. The symbol <j) denotes the char- 
acteristic vibrational temperature which is unique for each gas species. For N 2 we use 
(j) = 3395 0 K. For small temperatures we have C v ~ 5R/2 so that the vibrational degree 
of freedom only becomes excited when the temperature is on the order of magnitude of (f>. 
The coefficients of viscosity 77 and A = — 2r;/3 are determined from the Sutherland formula, 
reference [3] where 

<=10=T 3/2 

’ = -tT^T (10) 

where for iV 2 we use c\ = 1.488 * 2.16(10 -8 ), g c = 32.174, and c 2 = 184.0. The units of 
viscosity axe kg/m — sec when the temperature T is given in Rankine units. The quantity 
$ represents the dissipation function and is given by 


* = V(r 0 V y ) - V • V(r,j) = 

where are the viscous stress terms given by 

Ti j =r,(V i , j + V jti ) + \S ij V-V. 

The coefficient of thermal conductivity K is written, reference [4] 

ic = if rl + jr„ = ^ + ^ 


(ii) 


( 12 ) 


(13) 


where P r is the Prandtl number and S c is the Schmidt number. For N 2 at 0 0 K , we have 
the approximations P r ~ 0.71 and S c ~ 0.74. These values remain constant over a wide 
range of temperatures and produce the approximations 


K rt ~ 4 . 9377 R K v w l.35r]C vv . 


(14) 
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Computational Coordinates 


The basic equations (1) through (5) are written in the weak conservative form, 

Continuity 


dg 

dt 


1 d . i r \ 

+ -^ v ') + 


a( g v,) 

dz 


= 0 


( 15 ) 


Momentum 

I W) + +p- T rr]) + l z (ev r v, -r r ,)-^+ T f=0 (16) 

§(eK) + - r„)) + + P - r 2I ) = 0 (17) 


Energy 

np I n ^ 

-^- + -^(r[(^t+-P)V r -V r r r r -V 2 r„+g r ])+ — ((E t +P)V' 2 -V r r rz -V z T 22 +g z ) = 0 (18) 

Equation of State 


P = gRT (19) 

where E t is the total energy per unit volume 

E, = <?e + |(K 2 + V?) (20) 

and the stresses given by 


dV r 

r rr = 2^ + AV • V 

(21) 

dV z 

t zz = 2 V -^ + AV-V 

(22) 

r ee = 2rj— + XV-V 
r 

(23) 

(dV z dV r \ 
= r " = ’ ! (ar + dz) 

(24) 

v - v 

(25) 


and e is the internal energy per unit mass determined from the relation 


de = C v dT. 

The internal energy per unit mass e is given by equation (9). 


( 26 ) 
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The equations (15) through (19) axe to be solved over the solution domain 0 < z < 6, 
0 < r < f(z) where f(z) defines the shape of the nozzle. We introduce the dimensionless 
variables 

t 


r =8 


Z L 


* 0 * V 

Q — — V = — 


0o 


Vo 


t* = 


P* = 


L/V 0 

P 

QoVo 


K ~ 


Vr 


V 0 6/L 

T 

T* = — 
T 0 


v; = 


e = 


K 

V 0 

e 

W 


(27) 


where L,8 are characteristic lengths, Vo a characteristic velocity, po a characteristic density, 
r/o a characteristic viscosity, and To is a characteristic temperature. There then results the 
dimensionless equations 
Continuity 


v + 4£ (rVv:)+ S£S2. 0 


dt 


r* dr* 


dz* 


(28) 


Momentum 


J?(»*K*) + ^(r'WV; + - ^KJRe 1) 


<* 2 


4 - (p*v*v* - 

+ dz* r z 


L 2 


8 2 

L 2 P* 


L 2 


—t- /Re) - — — + - lil_ 

6 2 ^ r „ + 62 r , Re 


= 0 


Jrte'K*) + V.^(r*[e*v?v; - p ;,/Re)) 

+ jp(e' v : v : + P ‘~ = 0 


(29) 


(30) 


Energy 


^ + 4 ^(>-*[(-Er + OK* - V^/fie - T* t / Re + S J]) 


at* 


r * Q r * 

d 


+ + P‘)V; - V r * T* t j Re - V;r’JRe + = 0 


(31) 


Equation of State 


p* = e *PT’(r 0 /F 0 2 ) 


( 32 ) 


where E* is the scaled total energy per unit volume 



and the stresses given by 


dV* 

T *rr = V dr : + A*v* • R* 

(34) 

dV* 

r 2 * z = 2r 7 *^- + A*V*-y* 

(35) 

T* ee = 2r ] *^ r + \*V*-V* 
r * 

(36) 

- =T . =v -( av ' + ^ dv ‘) 
rz zr 1 V dr* L 2 dz* J 

(37) 

Id dV* 

r*dr* K r) dz* 

(38) 


with 

, L 2 -KT 0 dT* m -KT 0 dT* 

Qr ~ 6 2 ^ 0 V 0 2 Re dr* ?2 ~ rj 0 V 2 Re dz * 

where Re = QoVoL/r] 0 is the Reynolds number, A* = —2rj*/Z and E* = Et^PoV^). 
We write these dimensionless equations in the weak conservative form 


dU* 

dt* 


+ 


j__a_ 

r* dr* 


(r*G*) + 


dF* 

dz* 


+ —~H* = 0 


(39) 


where 



r p * i 

1 1 

Q'v; 1 

II 

* 

Q *v z * 

II 

* 

e'v;v; + £p> - g-KJRe 
q-v;v; - 

I 

L E* t J 

1 1 

(e; + p')v;- v;r; r /Re - £v;t ;jRe + J 


(40) 


e‘v; 


0 

e'v;v; - gr;,/Re 

H* = 

+ £ree/Re 

i"v;v; + P'- t ; jRe 

y 

0 

(E; + P')v; - V r 'r' r JRe - V^r'JRe + q • J 


0 


(41) 


The solution domain is now 0 < r* < f(Lz*)/8 and 0 < z* < b/L. 
The change of variable 


x = 



V = 


f(Lz*)/8 


converts the system of equations (39) to the form 


dU 

dt* 


+ 


dE 

dx 


dF 

+ ^- + H = 0 
dy 


(42) 


( 43 ) 
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V 


over the domain 0 < x < 1 and 0 < y < 1 where 


E = —F* F = - G* - Ly f F* 

b ' f(bx) U f(bx) * ’ 

H = H^-F* + + H*), U = U* 


f(bx) yf(bx ) 

Note that in the limit as y — ► 0 we have the result 


lim 

y—*° 


G* + H* 

y 


* dVT 

o 

Q K r dy 6 2 Re dy 


\(E* + p*)W Z _ Tj^d£_ _ Li J_ (V *dr^ . T * <>¥L) . M 
L y^t ^ r ) dy Re dy S 2 Re\ V z dy ^ T rz dy dy 

Using the vector of primitive variables 

v = Co\( 0 *,v;,v z *,t*) 

and the vector of computational variables 


( 44 ) 


(45) 


(46) 


v = Coi((>‘, e 'v;, e'v;,E;) (47) 

we obtain the first three primitive variables from the computational variables from the 
relations 


K = U x = Q* (48) 

v 2 = U 2 /U x = v; (49) 

V 3 = U 3 /U x =V z *. (50) 


The remaining primitive variable is determined using Newton’s iterative method on the 
system 

rp rp e (-^~n) gQ /cf\ 

n+1 “ Tn C v (T n ) (51) 

where e 0 = Vq(U 4 — ^V x (jjV^ + V 3 2 ))/Vi and e(T) is given by equation (9) with ^ = C v . 
After solving for T we calculate V 4 = T* = T /To- Conversely, we can construct the 
computational variables from the primitive variables using the relations 

U X =V X 

u 2 = v x *v 2 

U 3 = V X * F 3 (52) 

U 4 = V x * e(T 0 * Vi) + Vi * (^V 2 2 + V 3 2 )/2. 
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Operator Splitting 


The weak conservative form for the equations of motion in terms of the computational 
coordinates x, y are given by the system of equations (43). We can then define the operator 
L x as the numerical solution of the system 4^ + ^ = 0 given by 


Predictor: 

Corrector: 


V'ij =U'„ - £(K+,,i - Eh) 

V?j =\w,i + - f- x {EZ - EZ U )). 


Define the operator L y as the numerical solution of the system 4^ + ^ ~ 0 


Predictor: 

Corrector: 

Define the operator L as the numerical solution of the system 4^ + H = 0 as 

Predictor: £/** =U* , - AtH * , 

Corrector: V% =\ (V’j + U% - A tH%) 

where F*j = F(U?j), = H(Vf, J), etc. 

For the method of operator splitting we time march according to the sequence of 
operators 

v $ 2 = l,l > lll,l,ui 

and progress the solution from time nAt* to (n + 2)Af* where At* is selected to satisfy 
the Courant-Fredrich-Lewy CFL stability condition. The CFL stability condition is deter- 
mined by time steps A t z and A t r evaluated at all internal node points. These time steps 
are given by (reference [7]) 


Vh =U-J - - FTj) 

V"i =\(Vh + Wi **-.))• 


At, 
A t r 


lKi + v^r+i(^ + ^) 

f(bx)ki 

|7 r | + VlRT + l ^ f(bx)k x C p + ^TET 3 -) • 


From all such time steps we select the minimum time step At = {At 2 , At r } and then 

scale this real time to calculate the scaled time step At* = . 

Other methods to solve the system of equations (43) are Runge-Kutta methods and 
various implicit methods. Current research is trying to establish an efficient method for 
solving the system of equations (43) together with appropriate boundary conditions. 


7 



Two Temperature Model 

For our second model we assume a vibrational degree of freedom together with a 
combined rotational-translational degree of freedom. Each degree of freedom is assumed 
to follow a Boltzmann distribution and the energy content of each degree of freedom is 
characterized by temperatures T v and T respectively. As the gas passes through the nozzle 
there is a certain finite relaxation time r for the vibrational mode of excitation to achieve 
equilibrium with the rotational-translational mode of excitation. Define the quantities: 


rii Population density of ith energy level 
Energy per molecule of the ith level 
hi Time rate of change of n; due to collisions 
q t Heat flux 
q = qrr + q* 

where q rr is the heat flux due to rotational energy 

q* Heat flux due to energy excitation of all energy levels 
* Total energy from all energy levels 

<T = njejUj 

i ^ 

where U{ =Vi — V is the diffusion velocity of molecule in state i 
$ The Dissipation function. 

Let ge* = Yh n i e i denote the total energy per unit volume of the vibrational mode so 
that by integrating over the volume and surface of an arbitrary volume element we obtain 

J ge* dv = — J ^ n;e,Ui • dS + J ^ h,e, dv 


where dv is a volume element and dS is an area element of the control volume and ri,- are 
rate equations to be determined. Using the Gauss divergence theorem and interchanging 
the order of summation and integration there results 


d_ 

dt 


(ge*) + ]T VinMUi + V)) = ^ h t e t 
i i 

p\ 

+ Vfoe'V') = £ |A, e , - V(n, e ,(7.) 


( 53 ) 


Using the identities 


P + V ■ V( S e*) + V(ee'V) - S e‘V ■ V — V ■ V( S e") 
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together with the continuity equation and 


D De* Dg De* T -» 

—(ge ) = g— 1 - e — — = g— ge V • V 

Dt K e ’ *Dt Dt * Dt 


we write the vibrational energy equation as 

De * 


e^ = ^n iCi -v q 


Assuming linear harmonic oscillators and using the rate equations from Meador, et. al., 
reference [5], 

_^ n , e , = __ 


where the subscript e denotes equilibrium, the above equation becomes 

De* e* — e* - 

— — = — V • q*. 

Dt t g 

Letting * denote the vibrational mode then if T v exists, the vibrational energy equation 
can be represented as 


De v DT V 

= 6 


Dt 


Dt 


e„(T) - e v (T„) 




(54) 


The second energy equation is obtained by substituting e = e rt + e v into the energy 
equation (4) to obtain the form 


+ e v) + ■ V(e r * + c v ) + PV • V + V(g r ( + q v ) — $ -f gC vv X — gC vv X = 0 (55) 


where 


C VV X = 


e v {T) - e„(T„) 


(56) 


and then employing the vibrational equation (54) to obtain the coupled energy equations 


r\ 

Q ~W + gV ■ Ver< + PV • V + V • - $ + qC vv X =0 

r\ 

e-~ + t v • Ve„ + v • q, - nC m X =0 


(57) 


which reduce to the form 

DT -1 


Dt gCrt 
DT V . -1 


Dt gC vv 


(v • q rt — $ + PV • y + ecux) 
V • ft + 


( 58 ) 
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where 


r T r T / Ax'- 

S.(T)» / C ut ,dT = / i?( 

^ T v JT V \ / 


* v ' Jt v V Jt. \TJ (e*/T-iy 

The integral in equation (59) is used to calculate the quantity 


x = (e v (T) - e„(T„)). 

C/t)U T 


Integration produces the result 


T 2 1 - e~^/ Tv 

Y ± v c 

<j)T 1 — e~^/ T 


{exp [*(£-f)]-l}- 


The quantity gC vv X = j(e v (T) — e v (T v )) can be viewed as a coupling term for energy 
between the vibrational and rotational-translational modes. The other terms in the coupled 
equations (43) are given by 


C v — C vr t + C v v , C vr t — 5i?/2, C vv — R(<f)/T v ) 2 


(e*/ T » - l) 2 


e r , = -RT 


r Tv 

° v — I C v v (IT = 
J To 


R<f> R<j> 

e <t>/T v _ i e <t>/T 0 _ l 


x = r—MT) - e v (T v )) 

^vv T 

x . S (lz±XT\ U 


) { exp [* (s - f )] " *} 


-e-^yrTl^ T) J ‘J 

€t = -K rt VT (65) 

q v = -KvVT v (66) 

K ri = 4.93^i2 (67) 

K v = 1.3577C,,,, (68) 

c a T 3 / 2 

77 = — q- Sutherland’s formula (69) 

c x = (1.488) * 2.16(10 -8 ), g c = 32.174, c 2 = 184.0 (70) 

Here we have assumed that the coefficient for self diffusion between molecules in different 
internal states is a constant for all energy states. This in turn produces the above specific 
heats. 

The divergence of the heat flux terms are given by 


V • Qrt = ~K r 


'd 2 T ldT d 2 T\ dK rt fdT 
dr 2 + r dr + dz 2 J dT 

d 2 T v 1 dT v d 2 T v \ dK v \(dTl 


ldT v d 2 T v \ dK v (dT v y (dT v 
r dr + dz 2 ) dT v \ dr J + \dz 
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For the pressure we assume an equation of state for an ideal gas P = pRT. Following 
Meador et.al. [5] the relaxation time r for N 2 is taken as 


P(atm)T = 


3.2188(10 -12 ) 
I(T) sinh(<£/2T) 


(?) ex p(-t/ T > 


(73) 


where <j>, axe characteristic temperatures given by <f> = 3395 K, 6 = 3.2324(10 7 ) if, 
£ = 95.9 K , and 





TV 


1/2 


1 - 1/3 



(l — e 2 ,: ')exp(-(x + < + -C_)]<Ji 


(74) 


where 


c± 


'16 6Tx 
_ 27(f) 2 



(75) 


For the two temperature model we replace the single temperature equation by two 
temperature equations. We define the column vector of primitive variables 


V = col( 0 , V,, V„ T, T„) 


(76) 


and use the computational variables 


U = Col( g , <?V r , £rt, £.)• (77) 

This requires that we replace the temperature equations (58) by energy equations in con- 
servative form. The conservative form of the split energy equations are obtained by sub- 
stituting E t = E rt + E v into the equation (18) to obtain the equations 

+ + P)Vr - VrTrr - V z T rz + <j> rt ]) 

+ -^[(E rt + P)V Z — V r T rz — V z t zz + q zr t) + gC vv X = 0 
d -^ + ~ME.Vr + qrv)] + -^{EvVz + qzv] ~ QC VV X = 0 

where 

E„ = f?e rl E rl = j0i7r + !(V r 2 + K 2 ) 

Ev = Q6v E v = gR<f>/{e' l>/Tv - 1) 

These equations must be scaled and added to the continuity and momentum equations 
developed earlier. The resulting set of coupled equations are given by 


dU* 

dt* 


+ 


1 d 

r* dr* 


(r*G*) + 


dF* 

dz* 


+ —H* = 0 


l 
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where 


r p* i 


e'v 7 1 

e*v* 


e’v;v; + £p- - f? t-jr? 

q*v. : 

II 

# 

03 

e'v;v; - £t ;,/Re 

Kt 


(e;, + P‘)v; - v;r; r /Re - £v; T ;jR e + 

L E* J 


b;v; + ,;„ J 


F* = 


where 


# 

1 


0 

e'v;v; - jyr ;jRt 


—P* + ^ 2 Toe /Re 

o-v:v; + p- - r:jRt 

II 

* 

0 

(e;, + p-)v; - v;r;jRe - v;r ;,/Re + g ; r , 


r*Q* 

e;v; + q - v J 


—r*Q* 


Kt = 

Ert 

E* - Ev 

eoV?’ 

" SoV?’ 

* 

L q rt 

* Qzrt 

<lrt = 

6 eoV’ ’ 

q "‘ “ Sol? 

n* — 

E Qrv 

* Qzv 

<lrv — 

6 eoVf’ 

5 '” - SoV 0 3 


Q ~ V* ’ 

Cvu = C'««|t=t 0 t*, ** = *| 


t v =t 0 t; 

T = T 0 T* 


T, 

To 


The above equations involve the primitive variables 


V = Co\(Q,V r ,V„T,T v ) 


and conservative variables 


U = Col(^, gV r , gV z ,E rt , E v ). 
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Boundary Conditions 

The initial inlet temperature, inlet pressure, and inlet density, are given a priori. From 
these values we calculate the critical values V*,P*, g*,A*,T* from a one dimensional model 
for the flow and from these values determine the inlet velocity. The one dimensional nozzle 
values of pressure, temperature, density and velocity are calculated as a function of nozzle 
distance z. These values are then used as starting values for each model. We require a 
no slip boundary condition at the nozzle walls, symmetry with respect to the centerline, 
and normal derivative of the pressure to be zero at the walls. Extrapolated boundary 
conditions are then applied to the exit values and as well as extrapolated temperature 
values at the wall and centerline. The exit pressure is also initially assigned. 

Comparison with isentropic one-dimensional model 


For comparison purposes we also assume an isentropic process and calculate the results 
for a one- dimensional flow through the nozzle in the z direction where the area of the nozzle 
is a function of z. For an isentropic process we have 

dh - 

ds = — y*- = 0 , 


with h = e + RT , e = C V (T) dT and dh — C p dT. Consequently, 


, dT dP 
ds = C p — R~p~ = 0 


or 


Let 


[ p dP = [ T [Z+ (£\ 

Jpo p Jto 2 \t) 


2 P 4>/T 


( e */r_i)2 


dT 
T ' 


V = 


1 


e*/ T - 1 ’ 
and integrate by parts to obtain 


e^ T (j) 


HP I?.- 2 logT I?. + f e */T-i lr. + J Ta 


* 


T 2 e*! T - 1 


dT. 


( 78 ) 

( 79 ) 

( 80 ) 


( 81 ) 


In the last integral, let z = e^/ T — l,dz = -t^^-^dT so that 

r T 4 i r dz _ r dz r dz 

J To T*e*/T- l dl - J w z(z + 1) - J„ z + /„ z + 1 


and consequently 



±__i_ 

T2 e H T - l 


dT = log 


/ 1 - e-+/ T ° 

V 1 - 


( 82 ) 
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Therefore, we can calculate the pressure ratio as 


P (T \ 7/2 ( <f>/T t/T 0 \ 

P 0 ~\To) V 1 - ) 6XP \e*/ T - 1 e*/ T <> - 1 )' 

Since P = qRT we can write 

e f t/T <t,/T„ \ 

So J eXp \e*l T -1 e*m - l) ' 

where (f>/T 0 is treated as a parameter. 

In one dimension the energy equation can be written 

dh + V z dV z = 0 or C p dT + V z dV z = 0. 


Consequently, we may write 



V z dV 2 = - 



C p dT 


which integrates to 


V; - V; o = 7R(T 0 - T) + 2 Rj> 


t 4>/To _ 1 e 4>/T 


M' 


( 83 ) 


(84) 


(85) 


(86) 


By dividing by a 2 = 7 RT, the local speed of sound, the one dimensional mach number 
can be represented 



T ) 7 T 


Z<t>!To _ J e <t>/T 


b) 


(87) 


with M = V z /a. The mach number and one dimensional analysis is used to obtain an 
approximate solution to the more complicated two dimensional problem. Here 


a, 7 + 2wr)V/7y/7’ - 1 )- 2 

7 C„ -j + 2(J,IT)' 2 e‘ilT[ e ilT _ lr l' ' ^ 

The one dimensional continuity equation is given by 

AV z g = A*V; e * (89) 

where the * quantities represent those values at the throat of the nozzle where M — 1. 
That is, set M = 1 in equation (84), then solve the equations (84)(85) simultaneously for 
the value of <f>/T , treating (t>/T 0 as a parameter. This calculated value of <^>/T gives T = T* 
when M = 1 and consequently we can calculate the values of P*, g*, 7 *, V* = 7 *RT* at 
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this critical value of the temperature. The equation (86) can then be expressed in the 
following form involving the above critical parameters 


v: 

A_ _ V*Q * _ y/rlvF Rg*T* 
A* V z g v M __ RqT^t 

y/^RT y/yRT* 


A 1 ly*T* P* _ 1 frTP* 
a*~m\I 7 t pJf- ~ M]J 7 T* P 

A_ = l_ IrT^Po 
A* M Y 7 <f>T* P 0 P' 


(90) 


Knowing the critical values T*,<y*,P*,A* we can calculate the ratio T / <f> as a function of 
A/ A* which is a function of z, with To as a parameter. These one dimensional values are 
then used as starting values for the solution of the two dimensional non-isentropic nozzle 
problem. 
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APPENDIX A 
LIST OF SYMBOLS 


a 

A 
— ♦ 
b 

C v , Cvrt j Cvv 
Dij 

=! +i7 ' v 
e, e rt , e„ 

E rt ,E v ,E t 

Ti 


h 

k 

K, K rt , K v 
m = W/N a 
M 
Na 

9 ) Qrt > Qv 

P 


Q 

R 


r 



Speed of Sound [m/s] 

Cross sectional area [m 2 ] 

Body force per unit mass [Newton /Kg] 

Specific heat at constant volume [Joule/ Kg K] 
Rate of deformation tensor [s -1 ] 

Material or substantial derivative 

Energy per unit mass [Joule/Kg] 

Energy per unit volume [Joule /m 3 ] 

Planck’s constant [Joules] 

Enthalpy [Joule /m 3 ] 

Boltzmann’s constant [Joule/ K] 

Thermal conductivities \W/mK] 

Molecular mass [Kg] 

Mach number 

Avogadro’s number [mo/ -1 ] 

Energy flux [Joule/m 2 s] 

Pressure [Newton/m 2 ] 

External heat source per unit volume [Joule/m 3 ] 
Gas Constant [Joule/ Kg K] 

Radial distance [m] 

Entropy per unit volume [Joule/m 3 K] 

Time [s] 

Temperatures [K] 



V 

V r ,V z 

w 

X 

x,y 

z 

V 
A 

Q 

T 

T a 


V 



$ 


Velocity [m / s] 

Velocity components [m/s] 

Molecular weight of iV 2 [Kg /mol] 
Coupling term [K/s] 

Computational coordinates 
axial distance [m] 

Viscosity coefficient [Kg /ms] 

Second coefficient of viscosity [Kg /ms] 
Density [Kg/m 3 ] 

Relaxation time [s] 

Stress tensor [Newton/m 2 ] 
Characteristic temperatures [K] 
Frequency [s -1 ] 

Ratio of specific heats 

Dissipation function [ Joule/m 3 s] 
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